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DPI 7066 

DEMANDE DE BREVET ^INVENTION 
METHODE DE CALAGE D'UN PUITS DE FORAGE DANS UN BLOC SISIWIQUE 



Invention de : Naamen KESKES 

Frederic MANSANNE 



SOCIETE ANONYME D1TE : 
TOTALFINAELF 

ABREGE DESCRIPTIF 

La m6thode a pour objet de determiner la meilleure correlation possible pour 
des enregistrements 6lectriques acquis dans un puits de forage avec des 
enregistrements sismiques acquis dans un volume de sous-sol contenant le puits. Le 
but est d'obtenir le meilleur calage possible du puits de forage dans le bloc sismique. 

Pour ce faire, la methode definit dans un voisinage de la trajectoire du puits 
un ensemble de reseaux en couches dont les entries sont les signaux sismiqu s et 
les sorties sont les signaux 6lectriques et elle attribue une mesure de quality a leur 
convergence. 

La meilleure qualite de convergence permet alors de determiner le meilleur 
lieu dans le voisinage qui produit le meilleur calage. 
(Figure 3 pourabrege) 
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DESCRIPTION 

La presente invention concerne une methode de calage d'un puits de forage dans un 
bloc sismique multi-attributs, e'est-a-dire qu'elle vise a obtenir la meilleure position 
d'une trajectoire de forage dans le bloc sismique. 

5 Etat de Tart 

La mdthode selon I'invention, s'applique sur des blocs sismiques. Un bloc sismique 
est forme par la juxtaposition de traces sismiques dans un volume. On appelle 
"coupe sismique" une tranche de bloc sismique. Une coupe sismique offre une vue 
de la juxtaposition des traces sismiques contenues dans le plan de coupe. Dans une 
10 image sismique, I'intensite lumineuse d'un pixel est proportionnelle a I'attribut 
sismique represents. 

Le calage d'un puits de forage dans un domaine sismique a deja fait I'objet de 
quelques travaux. 

En particulier, le brevet FR-A-2 737 909 propose une methode pour positionner une 
15 trajectoire de forage sur une section sismique par moderation dans un domaine 
somme synthetique des lieux de la trajectoire pour lesquels on a connaissance du 
couple (pendage, azimut) puis par migration en temps de I'image synthetique. 

Cette methode permet de placer le puits sur une coupe sismique par une methode 
de calcul deterministe. 

20 D'autres methodes de calcul ont ete proposees. Cependant, toutes sont 
deterministes. 

Apport de I'invention 

L'approche de la methode selon I'invention est de nature non-deterministe. On admet 
que toute mesure est entachee d'erreur et on s'attache a determiner le meilleur 
25 calage possible par un balayage de toutes les possibilites de calage pour en 
selectionner la meilleure. 

La methode selon I'invention determine la valeur du calage par la quality de 
convergence en phase d'apprentissage d'un ensemble de reseaux de neurones. 
La nouveaute reside dans le fait que les reseaux de neurones sont non seulement 
30 utilises comme estimateurs universels mais aussi comme instruments de mesure. 
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Ainsi chaque r6seau produit une mesure de "distance" de calage et le meilleur 
calage est determine par la distance minimum. 

Cette operation est analogue a la selection de la distance minimum dans un 
ensemble de mesures de distance prises avec un quelconque instrument de mesure, 
5 tel qu'une regie graduee. 

Definition de ['invention 

La pr6sente invention a pour objet une methode de calage, dans un espace 
presentant deux dimensions spatiales et une dimension temporeile, d'un signal S 
mesure en des positions U affect6es d'une incertitude, sur un ensemble de N 
10 signaux mesur6s en des positions determinees, les N + 1 signaux ayant leur origine 
temporeile dans un meme plan, ladite methode consistant a : 

- r6-6chantillonner les N + 1 signaux afin de les placer tous dans une plage 
d'echantillonnage identique, 

- filtrer le signal S afin de le placer dans une plage de frequences identique a celle 
15 des N signaux, 

et caracterisee en ce qu'elle consiste en outre a : 

- definir pour chaque position U associee aux mesures du signal S un m6me 
voisinage de Heux V dans I'espace spatio-temporel centre sur la position U, 

- fabriquer un r6seau de neurones en couches RN V pour chaque lieu V dans le 
20 voisinage de U, chaque r6seau ayant un vecteur d'entree de dimension N associe 

aux mesures des N signaux et une sortie scalaire associee d une mesure du signal 
S, 

- pour chaque r£seau de neurones RN\ definir un ensemble d'apprentissage tel que 
les entrees sont la collection de tous les vecteurs de mesures des N signaux situ6s 

25 aux lieux V et les sorties sont la collection des valeurs du signal S aux positions U 
pourtoutes les positions U, 

- fixer un nombre derations predetermine Nit pour tous les reseaux de neurones et 
lancer les phases d'apprentissage de tous les reseaux, 
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- pour chaque reseau de neurones RN V , calculer la valeur de i'integrale z v de la 
fonction donnant I'erreur commise par le reseau a chaque iteration, de I'iteration 1 a 
I'iteration Nit, 

- pour chaque position spatiale V k en surface du voisinage de coordonnees (Xk.yx .to), 
5 selectionner dans la dimension temps le couple de lieux (V1 k (x k ,y kl ti), V2 k (x k ,y k ,t 2 )) 

du voisinage qui correspondent aux deux plus petits minima locaux des deux 
integrates (T, v \ 7 Z v2 k ), 

- pour chaque position spatiale V k en surface du voisinage, retenir parmi les deux 
positions V1 k (x k ,y kl t 1 ) et V2 k (x k ,y k ,t 2 ) la position V m , pour laquelle le signal estime par 

10 les reseaux de neurones respectifs RN v1 k et RN V2 k presente une variance maximale, 

- choisir parmi les V m la position pour laquelle I'integrale E v m est minimum. 
Figures 

La figure 1 represente un signal electrique positionne dans un bloc sismique. 

La figure 2 illustre le resultat du filtrage d'un signal electrique qui ajuste sa plage de 
15 frequences a celle des signaux sismiques. 

La figure 3 represente un schema de principe de fonctionnement des reseaux de 
neurones dans la methode selon ('invention. 

La figure 4 montre quatre fonctions d'erreurs associees aux phases d'apprentissage 
de quatre reseaux en couches. 

20 La figure 5 illustre le principe de selection du meilleur decalage temporel. 

La figure 6 presente un reseau de neurones dans lequel les entrees sont issues 
d'une fenetre temporelle prise sur les signaux sismiques. 

La figure 7 montre le resultat de I'estimation d'un enregistrement electrique realisee 
par le reseau de neurones et appliquee au meilleur point de calage. 

25 Description complete 

La methode selon I'invention est une methode de calage d'un puits de forage dans 
un bloc sismique. En se referant aux figures, un mode de mise en oeuvre de cette 
methode est donne ci-apres. 
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Une acquisition sismique permet de fabriquer un bloc sismique dont le but est de 
representor le volume de sous-sol parcouru par Tonde sismique. Un bloc sismique 
est constitu6 d f un ensemble de traces sismiques somm6es, c'est-&-dire pretraitees 
pour obtenir la meilleure vision possible des reflecteurs sous-terrains. Les r6flecteurs 
5 sont des interfaces entre deux milieux presentant des vitesses differentes pour I'onde 
sismique. Les reflecteurs represented generalement des limites de couches 
geologiques. 

Soit un bloc sismique B. Le bloc sismique est defini par une matrice M & trois 
dimensions. Deux dimensions du bloc B sont H6es a deux axes spatiaux horizontaux, 
10 appeles commun6ment (Ox, Oy). La troisieme dimension du bloc B est Ii6e a un axe 
temporel vertical, appele commun6ment Ot. 

Le bloc sismique 6tant une image du sous-sol, chaque element de la matrice M est 
assocte a un point de coordonnees (xbi, yb lt tbj) dans le sous-sol. La composante 
temps tbi est Ii6e a la profondeur zbj par une relation du type : 

15 zb» = v(tb,) * to* (1) 

ou v(tb|) est une fonction qui donne la vitesse de Tonde sismique en fonction de la 
profondeur. 

Chaque element de la matrice M contient un vecteur d r attributs sismiques (An, A^,... 
A io ) au point de coordonnees (xbi, yb jt tbO- Par exemple, les attributs (A*, A i2f ... A } „) 
20 mesurent I'^nergie de I'onde sismique, aussi appelee amplitude sismique. 

On appelle deport la distance qui separe le couple Gmetteur-recepteur associ6 a une 
trace sismique. Une mesure d'amplitude sismique est alors determine par la 
sommation des mesures d'amplitudes sismiques selon des classes de deports 
predetermines. Typiquement, on sommera les traces selon une partition en trois 
25 classes de deports : deports proches, d6ports lointains et deports tr6s lointains. 

Soit un puits de forage F for6 dans le volume de sous-sol couvert par le bloc 
sismique B. Le puits F parcourt le sous-sol selon une trajectoire de forage definie par 
un ensemble de points de coordonn6es spatiales (xfi,yf|,zf|). 

Tout forage de puits d'exploration p6troltere donne lieu & I'enregistrement de 
30 nombreuses mesures 6lectriques dans la profondeur du purts. Pour chaque point de 
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coordonn&es (xfi,yfj,zfj) de la trajectoire du puits de forage, on acquiert ainsi un 
ensemble de mesures electriques (E-| F E2,-. E n ). 

La dimension des profondeurs zf, est I'unite de longueur ; les mesures sont prises en 
metres ou en pieds. Une premiere operation pour placer le puits dans le bloc 
5 sismique est de transformer la grandeur longueur en grandeur temps, Cette 
transformation s'opdre aisement grace 3 Ja relation (1) puisque Ton peut mesurer les 
vitesses de propagation de Fonde sismique dans le sous-sol au voisinage du puits, 
ces vitesses etant directement mesurees dans le puits de forage. 

Le resultat de la transformation est une trajectoire de puits temporelle definie par un 
10 ensemble de points U de coordonnees spatiales (xfi,yfi f tfj). 

On accroche ensuite le puits dans le bloc sismique par les coordonnees de depart de 
sa trajectoire ou coordonnees de surface (xfo, yfo, tf 0 ). 

La figure 1 schematise une portion d'un bloc sismique (cube sismique) presentant 
des pas en espace dx et dy dans les deux dimensions horizontals et un pas en 
15 temps dt dans la dimension verticale. Un signal electrique (Log) est place dans la 
portion du bloc. 

Comme les enregistrements electriques dans les puits de forage sont generalement 
effectues tous les pieds et que la resolution d'un signal sismique, fonction de sa 
bande passante, est de I'ordre de quelques metres, il convient dans un premier 

20 temps de re-6chantillonner les signaux electriques et sismiques de sorte que leur pas 
d'echantillonnage dans la dimension temps soient similaires. Pour ce faire, on r6- 
6chantillonne Tensemble des mesures electriques du puits F et on sur-echantillonne 
Tensemble des mesures d'amplitude sismique du bloc B. Dans un deuxieme temps, 
on applique un filtre sur les signaux electriques pour les amener dans une plage de 

25 frequences identique £ celle des signaux sismiques. 

La figure 2 montre un exemple de filtrage d'une mesure electrique. Le rectangle de 
gauche intitule "Log" represente Tamplitude d'un signal electrique 6chantillonne en 
temps apres application de la transformation (1). Le rectangle de droite intitule "Log 
Filtre" represente le signal electrique dont les frequences ont ete amenees dans le 
30 domaine des frequences sismiques par application de filtres passe-haut et passe- 
bas. Le rectangle superieur intitule "Spectre en Frequence", dont I'axe horizontal 
gradue des frequences et I'axe vertical des amplitudes, montre les deux spectres de 
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frequences associ£s £ une trace sismique (Near) et au signal electrique filtr6 (Log 
filtre). On peut noter que le signal electrique qui contenait 3 Torigine beaucoup plus 
de hautes frequences que la trace sismique se trouve maintenant ramene dans (e 
domaine frequentiel du signal sismique. 

5 On definit alors pour chaque point U constituant la trajectoire du puits et assocte h 
une mesure electrique Ei un voisinage spatio-temporel dans le bloc sismique B sous 
la forme d'un parallelepip&de centre sur le point U. A I'instar du bloc sismique B, le 
parallel£pip£de possede deux faces spatiales horizontals et une dimension 
temporelle verticale, orthogonale a la face spatiale. Soit Dvx le nombre de lieux du 
10 voisinage sur Paxe Ox f Dvy le nombre de lieux du voisinage sur l'axe Oy et Dvt le 
nombre de lieux du voisinage sur I'axe Ot. 

On choisira, par exemple, une dimension du parallelepipede de 5 6chantillons de 
part et d'autre du point U dans chaque dimension. Dvx, Dvy et Dvt prennent, dans ce 
cas, tous pour valeur 11 et le voisinage est d£fini comme un cube contenant 1331 
15 (11x11x11) points de mesures sismiques. Le meme voisinage s'applique a tous 
les points de la trajectoire, c'est-a-dire par extension h la trajectoire elle-mSme. Le 
voisinage est defini par un ensemble de points V. On appelle Npv le nombre de 
points inclus dans le voisinage. Npv est 6gal au produit Dvx x Dvy x Dvt et vaut 1331 
dans rexernple precedent 

20 On fabrique ensuite Npv reseaux de neurones en couches. Pour un cube de 1331 
points, on fabrique ainsi 1331 reseaux de neurones. Le vecteur d'entree de chaque 
reseau est un vecteur de mesures d'amplitude sismique dont chaque composante 
appartient a une classe de deports differente. La sortie du reseau se reduit £ un 
scalaire representant une mesure Electrique E». 

25 La figure 3 offre un schema de principe. Dans cet exemple, les reseaux de neurones 
possedent deux entrees et une sortie. Le rectangle de gauche intitule "Near" montre 
les valeurs prises par la premiere entree du reseau de neurones dans I'ensemble 
d'apprentissage. Le rectangle central intitule "Far" montre les valeurs prises par la 
seconde entree du reseau de neurones dans I'ensemble d'apprentissage. Le 

30 rectangle de droite intitule "Log" montre les valeurs prises par la sortie du reseau de 
neurones dans I'ensemble d'apprentissage. Un exemple est done constitu6 par deux 
valeurs d'entr6e prises a un temps (Time) determine sur les signaux "Near" t "Far" 
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et une valeur de sortie prise au meme temps sur le signal "Log". L'ellipse intitul6e 
"RN" est une representation symbolique du reseau de neurones entre entrees et 
sortie. 

Les r^seaux de neurones ont alors pour but d'apprendre une relation entre la mesure 
5 electrique E" et les mesures sismiques (Ai\ A 2 V # ... A n v ), c'est-a-dlre une fonction f 
telle que : 

E^ = f(A 1 v f A 2 V ,... A„ v ) (2) 

Une mesure indicee par u indique que sa valeur est prise sur la trajectoire du puits et 
une mesure indic6e par v indique que sa valeur est prise dans le voisinage. 

10 Les reseaux de neurones pr6sentent tous la meme architecture, c'est-a-dire le meme 
nombre de couches cach6es, le m§me nombre de neurones par couche, le meme 
nombre de connections et les memes conditions initiates. 

On choisit une position V dans le voisinage. La mise en aeuvre du r6seau de 
neurones RN V assocte & la position V s'effectue de la fagon suivante. Pour tous les 
15 points U de la trajectoire du puits, on fabrique un ensemble d'apprentissage pour le 
reseau RN V en collectant tous les vecteurs de mesures d'amplitude sismique situes a 
la position V. L'ensemble d'apprentissage pour le reseau de neurones RN V est ainsi 
constitue d'autant d'exemples qu'il y a de points U sur la trajectoire. 

On choisit un nombre d'iterations pr6d6termfne Nit et on lance la phase 
20 d'apprentissage du r6seau de neurones RN V par r6tro-propagation du gradient. 

Lorsque les Nit iterations sont effectu6es, on calcule la valeur de Pint6grale E v de la 
fonction donnant I'erreur quadratique commise par le reseau 5 chaque it6ration, de 
('iteration 1 3 Titration Nit. 

i 

On choisit ensuite une nouvelle position V dans le voisinage de la trajectoire et on 
25 accomplit les etapes precedentes avec le r§seau de neurones RN V ' associe a la 
position V pour obtenir une nouvelle integrate £ v '. 

On renouvelle ainsi Toperation jusqu'& ce que tous les points du voisinage aient ete 
balay6s. Pour un cube de 1331 points, on obtient ainsi 1331 valeurs d'integrale. 

En machine, ces operations peuvent se r6aliser de fagon sequentielle avec le m§me 
30 reseau de neurones reinitialise 5 chaque changement de point dans le voisinage. 
Cependant, il est recommande de mettre en oeuvre cette technique en parallele dans 
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des machines & multi-processeurs organises en fermes de calcul en fabriquant 
autant de reseaux de neurones que la capacite des configurations-machines le 
permet L'id6al serait de pouvoir obtenir autant de reseaux de neurones que de 
points dans le voisinage. 

5 Tous les points du voisinage ayant 6t6 parcourus, pour chaque position spatiale V k 
en surface du voisinage de coordonndes (Xk f y k to) au nombre de Dvx x Dvy, on retient 
dans la dimension temps le couple de lieux (V1 k (x kl y kl ti),V2k(Xktykft2)) du voisinage 
qui correspondent aux deux plus petits minima locaux des deux integrates (Z v1 k , 2^). 

La figure 4 montre quatre fonctions d'erreurs associees aux phases d'apprentissage 
0 de quatre reseaux de neurones identiques. Les axes horizontaux sont gradu6s en 
iterations et les axes verticaux donnent les erreurs quadratiques commises par les 
reseaux lors de leur estimation du signal electrique. Les valeurs des quatre integrates 
associees aux r&seaux sont egales aux surfaces hachurees. Dans cet exemple, ce 
sont les reseaux identifies par les triplets (dx k , dy kf dt|) et (dxi<, dy k , dt 2 ) qui 
5 presentent les deux plus petits minima locaux des integrates. 

On calcuie ensuite la variance du signal estim6 par les reseaux de neurones RN Vi k et 
RN^k et on retient dans le couple (V1 k (x k> yk l ti) f V2k(x k ,yk.t 2 )) la position V m pour 
laquelle la variance est maximale. On d£finit ainsi un ensemble compost de Dvx x 
Dvy lieux V m dans le voisinage. 

La figure 5 indique comment se fait la selection de la position V m a partir d s 
positions V k . Le rectangle sup6rieur donne les valeurs d'intdgrales en fonction du 
deplacement sur Faxe temporel pour un V k choisi. Dans cet exemple, le voisinage 
tempore! est fix6 & 20 6chantillons temporels de part et d f autre d'un point de la 
trajectoire. Les decalages DE1 et DE2 correspondent aux deux plus petits minima 
locaux des integrates. Le rectangle inferieur represente les valeurs des variances des 
deux signaux estimes par les r£seaux de neurones assoctes aux decalages DE1 et 
DE2. On retient le d£calage DE1 comma meilleur calage V m pour la position de 
surface V k . 

Enfin, on choisit parmi tous les Dvx x Dvy lieux V m , la position V ca i pour laquelle z v m 
est minimum. La position V^i est la position du voisinage qui assure le meilleur 
calage du signal electrique sur les signaux sismiques. 
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Un autre mode de realisation consiste £ ne pas renseigner les reseaux de neurones 
par une mesure unique par classe de deports mais par un ensemble de M valeurs 
prises dans une fenetre d'6chantillonnage on temps sur la trace, laquelle fenetre est 
centrde sur un point V du voisinage. On prendra, par exemple, une fenetre 
temporelle 6gale a 8 echantillons de part et d'autre du point V, c'est-^-dire 17 valeurs 
de temps. 

La figure 6 pr6sente un tel reseau pour fequel la fenetre temporelle choisie est de 2 
echantillons de part et d'autre du point V pour des raisons de simplicity. Le reseau 
possede de ce fait 5 entries par classe de deports. La figure montre deux classes de 
deports proche (Near) et lointain (Far) qui gen6rent 2x5 entrees. 

Ceci modifie I'architecture du r6seau qui comporte alors un vecteur d'entree dont la 
dimension passe deNaNxM. La sortie du r6seau demeure scalaire. 

La m6thode selon Pinvention s'applique ensuite de fagon identique et conduit a la 
selection du meilleur point du voisinage comme meilleur calage possible. 

Cette technique de fenetrage a I'avantage de permettre au reseau d f effectuer une 
d6convolution sur la base de la dur&e de I'impulsion sismique 6mise dans le sous- 
sol. La qualite des r6sultats s'en trouve alors grandement amelior6e. 

La figure 7 montre le resultat de ^estimation de I'enregistrement efectrique 
(Prediction) realisee par le r§seau de neurones a partir des signaux sismiques et 
appliqu6e au meilleur point de calage. Le trace en pointille (Puits) represente le 
signal electrique effectivement mesur6. On notera la tres grande qualite de 
(•estimation qui atteste de la qualite du calage d6termin6 par la methode selon 
('invention. 
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REVENDICATIONS 

1. M6thode de calage, dans un espace pr£sentant deux dimensions spatiales et 
une dimension temporelle, d'un signal S mesur£ en des positions U affectees 
d'une incertitude, sur un ensemble de N signaux mesur6s en des positions 
5 d6terminees, les N + 1 signaux ayant leur origine temporelle dans un m§me plan, 

ladite methode consistant a : 

- r6-6chantillonner les N + 1 signaux afin de les placer tous dans une plage 
d'6chantillonnage identique, 

- filtrer le signal S afin de le placer dans une plage de frequences identique a 
10 celle des N signaux, 

et caracterisee en ce qu'elle consiste en outre £ : 

- definir pour chaque position U associee aux mesures du signal S un meme 
voisinage de lieux V dans Tespace spatio-temporel centre sur la position U f 

- fabriquer un r6seau de neurones en couches RN V pour chaque lieu V dans le 
15 voisinage de U, chaque reseau ayant un vecteur d'entr6e de dimension N 

associe aux mesures des N signaux et une sortie scalaire associee a une 
mesure du signal S, 

- pour chaque r6seau de neurones RN\ d6finir un ensemble d'apprentissage tel 
que les entrees sont la collection de tous les vecteurs de mesures des N signaux 

20 situ§s aux lieux V et les sorties sont la collection des valeurs du signal S aux 

positions U pourtoutes les positions U, 

-fixer un nombre derations predetermine Nit pour tous les reseaux de neurones 
et lancer les phases d'apprentissage de tous les r6seaux, 

- pour chaque reseau de neurones RN V , calculer la valeur de Tint6grale S v de la 
25 fonction donnant I'erreur comrnise par le r6seau a chaque iteration, de Titeration 

1 a ['iteration Nit, 

- pour chaque position spatiale V k en surface du voisinage de coordonnees 
(Xk.ykto), selectionner dans la dimension temps le couple de lieux (V1k(Xk,yk,ti), 
V2 k (Xk,yk.t 2 )) du voisinage qui correspondent aux deux plus petits minima locaux 

30 des deux integrates (E v \, £ v2 k), 
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- pour chaque position spatiale V k en surface du voisinage, retentr parmi les deux 
positions V1 k (x k ,yk,ti) et V2k(x k ,y kt t 2 ) la position V m , pour laquelle le signal estime 
par les reseaux de neurones respectifs RN V \ et RN V2 k presente une variance 
maximale, 

- choisir parmi les V m la position V ca i pour laquelle integrate z v m est minimum, 

2. IVtethode selon la revendication 1, caracterisee en ce que la mise en aeuvre des 
reseaux de neurones consiste a : 

- definir pour chaque position U assoctee aux mesures du signal S un meme 
voisinage de lieux V dans I'espace spatio-temporel centr6 sur la position U f 

- fabriquer un reseau de neurones en couches pour chaque lieu V dans le 
voisinage de U, chaque reseau ayant un vecteur d'entree de dimension N x M 
associe aux mesures sur une fenetre temporelle de taille M centree en V des N 
signaux et une sortie scalaire associee a une valeur du signal S f 

- pour chaque reseau de neurones, d6finir un ensemble d'apprentissage tel que 
les entries sont la collection de tous les vecteurs de mesures prises dans une 
fenetre temporelle de taille M centrSe en V pour les N signaux et les sorties sont 
la collection des valeurs du signal S aux positions U pour toutes les positions U, 

- fixer un nombre d'iterations pred§termin6 Nit pour tous les reseaux de neurones 
et lancer les phases d'apprentissage de tous les reseaux, 

- pour chaque reseau de neurones RN V , calculer la valeur de I'integrale S v de la 
fonction donnant I'erreur commise par le r6seau a chaque iteration, de Hteration 
1 a ^iteration Nit, 

- pour chaque position spatiale V k en surface du voisinage de coordonnees 
(x k ,y k t 0 ), s6lectionner dans la dimension temps le couple de lieux (V1 k (x k> y kf t|), 
V2 k (x«,y k ,t 2 )) du voisinage qui correspondent aux deux plus petits minima locaux 
des deux integrates (£ v1 ki E^O, 

- pour chaque position spatiale V k en surface du voisinage, retenir parmi les deux 
positions V1 k (x k ,y k ,t|) et V2 k (Xk ( yk,t 2 ) la position V m , pour laquelle le signal estim6 
par les reseaux de neurones respectifs RN v1 k et RN V2 k presente une variance 
maximale, 
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- choisir parmi les V m la position V ca) pour laquelle I'integrale 2 v m est minimum. 
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